Spin- Wave Description of Nuclear Spin-Lattice Relaxation in Mn^O^ Acetate 
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In response to recent nuclear-magnetic-resonance (NMR) measurements on the molecular cluster 
Mni20i2 acetate, we study the nuclear spin-lattice relaxation rate 1/Ti developing a modified spin- 
wave theory. Our microscopic new approach, which is distinct from previous macroscopic treatments 
of the cluster as a rigid spin of S = 10, not only excellently interprets the observed temperature and 
applied-field dependences of 1/Ti for 55 Mn nuclei but also strongly supports the 13 C NMR evidence 
for spin derealization over the entire molecule. 

PACS numbers: 76.50.+g, 05.30.Jp, 75.50.Xx 
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Mesoscopic magnetism jy is one of the hot topics in 
materials science, where we can observe a quantum-to- 
classical crossover on the way from molecular to bulk 
magnets. Metal-ion magnetic clusters are thus interest- 
ing and among others is [Mni20i2(CH 3 COO)i6(H 2 0)4] 
H (hereafter abbreviated as Mnl2), for which quantum 
tunneling of the magnetization ^|JJ] was observed for the 
first time. There are three symmetry-inequivalent Mn 
sites in the Mnl2 cluster (see Fig. m. The four inner 
Mn 4+ spins and the eight outer Mn 3 " spins are directed 
antiparallel to each other and exhibit a novel ground state 
of total spin S = 10 ||. Resonant magnetization tunnel- 
ing in such high-spin molecules can be an evidence for 
the validity of quantum mechanical approaches at the 
nanometer scale. 

The simplest Hamiltonian for the Mnl2 cluster in a 
field may be given by H = -D(S Z ) 2 - gfi B S ■ H , 
where the molecular cluster is strictly treated as a rigid 
S = 10 object with single-axis magnetic anisotropy. 
While such a macroscopic treatment of the molecule 
interprets well quantum relaxation of the magnetiza- 
tion gj^] , recent electron-paramagnetic-resonance (EPR) 
measurements || suggest a possible breakdown of the 
spin-10 description. Although a microscopic treatment 
of each Mn moment is necessary for further under- 
standing of nanoscale magnets and is interesting in it- 
self, the total spin states in the Mnl2 cluster is too 
large even for modern computers to directly handle. 
Thus it is an idea J§ that the Mn 3+ -Mn 4+ pairs con- 
nected by the strongest exchange interaction J\ con- 
struct composite spins \. This eight-spin scheme ex- 
plains well the magnetization JBlJlG] , inelastic-neutron- 
scattering 0Jl^l , and EPR |p^ , |l3| measurements at suf- 
ficiently low temperatures. However, a recent skillful 
numerical-diagonalization study has shown that even 
the ground state is quite sensitive to the less dominant 
exchange interactions J2, J3, and J4 as well, throwing 
doubt on a parameter assignment within the eight-spin 
scheme |J, J\ ~ —150cm -1 , J2 ~ J3 ~ —60cm" 1 , and 
J 4 ~ 0. 

It is also unfortunate for this eight-spin model to be less 
applicable to nuclear-magnetic-resonance (NMR) mea- 
surements, which have vigorously been carried out for the 



Mnl2 cluster fT5]-|20f in recent years. Although the nu- 
clear spin-lattice relaxation time T\ can serve as a probe 
to the electron spin dynamics within the molecule, it is 
not yet interpreted beyond the spin-10 description. The 
early NMR studies on the Mnl2 cluster were performed 
by using proton [jl5|,0 or deuteron jl(| nuclei as probes 
and therefore the interpretation was plagued by the aver- 
aging effect over the numerous protons and the weakness 
of the hyperfine coupling to the Mn moments. In or- 
der to obtain more direct information on the local mag- 
netic properties, 55 Mn NMR measurements jl^jljj have 
recently been performed. There has also appeared an 
elaborate 13 C NMR evidence of the paramagnetic 
spin density of the cluster being delocalized over the en- 
tire molecule. Now we make our first attempt to interpret 
microscopically the nuclear spin-lattice relaxation-time 
measurements, developing a modified spin-wave theory 
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We introduce a microscopic Hamiltonian for the Mnl2 
cluster as 
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Y^Ji^i ■ Si + 2J 2 { Sl ■ Si + Si ■ si +1 ) 

1=1 



-2J 3 (si ■ s [+1 + -si ■ s l+2 ) + 2J A {Si ■ Si + S r S l+ x) 



+D2(Sn 2 



D 3 (Sty+9»BH(sl + Sf + S?) 



(i) 



© Mn(l): Mn 4+ (s=3/2) 
• Mn(2): Mn 3+ (S=2) 
O Mn(3): Mn 3+ (5=2) 
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FIG. 1. Schematic plot of the Mnl2 cluster. Inequivalent 
sites Mn(l), Mn(2), and Mn(3) are occupied by Mn 4+ ions 
with 8=5, Mn 3+ ions with 5 = 2, and Mn 3+ ions with 
S — 2, respectively. There exist four types of exchange inter- 
action between them: Ji, J2, J3, and J4, which are drawn by 
thick solid, dashed, thin solid, and dotted lines, respectively. 
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where s/, Si, and Si are the spin operators for the 
Mn(l) (spin § = s), Mn(2) (spin 2 = S), and 
Mn(3) (spin 2 = S) sites in the Ith unit, and N, the 
number of the elementary units, is equal to 4. The 
strongest exchange interaction J\ and the next lead- 
ing antiferromagnetic interaction J 2 were reliably esti- 
mated at — 150 cm -1 and — 60 cm -1 , respectively, but 
the rest of the parameters remain to be fixed. Taking 
account of previous investigations, we adopt two con- 
trasting parameter sets: (a) J\ = — 150 cm -1 , J 2 = 
-60cm- 1 , J 3 = 60cm- 1 , J 4 = 30cm- 1 ppl; (b) 



J\ = —150cm- 1 , J 2 = —60cm- 1 , J 3 = — 30cm , J4 = 
30 cm- 1 PJI^|. As for the anisotropy parameters, there is 
much less information. When the molecule is treated as a 
rigid spin-10 object, the macroscopic uniaxial crystalline 
anisotropy parameter D is determined so as to fit the 
zero-field separation between the M — ±10 and M = ±9 
levels, which is about 14 K |&p|. Hence it is natural to 
choose the local single-ion anisotropy parameters, which 
describe the J aim- Teller-distorted Mn 3+ ions |23ll , within 



the same scheme 12 . Setting D 2 and D 3 both equal to 
1.5 cm -1 , we obtain the excitation energy of 13.7 K in the 
following spin- wave treatment. The g factors are all set 
equal to 2 

We consider a spin-wave treatment of the Hamilto- 
nian (Q) introducing the bosonic operators for the spin 
deviation in each sublattice via sf = —s + al-.ai^, 
,s ; + = \/2sa[ 1 ; Sf — S - a\ 2 a^ 2 , S* = V2Sai, 2 ; 

Sf = S — a l3 ai,3, 5 ; + = \/2Sai, 3 . We carry out the 
Bogoliubov transformation in the momentum space, 



&k,2 
1k,3 



-^22 (k) 6 fc ,2 
-^32 (k) b k ,2 



-tpi 3 (k) b K3 , 

-r 23 (k)bi 3 , 
-r 33 (k)bi 3 , 



so as to reach the diagonal Hamiltonian 
E E u 3(k)b\ d b k 



H = E S 
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(2) 



(3) 



k 3=1,2,3 



where -B g = 8Ss(,h + 2J 2 ) - 12s 2 J 3 - 16S 2 J±- 4S 2 (D 2 + 
D 3 ) — g^H(2S — s). The numerically calculated disper- 
sion relations u>i{k) for the two parameter sets are shown 
in Fig. 0. The lowest-lying ferromagnetic (o) and the 
antiferromagnetic (x) branches are both sensitive to J 3 , 
while the second ferromagnetic (o) branch exhibits little 
dependence on J 3 . 

The core idea of the so-called modified spin-wave the- 
ory is summarized as constructing reliable thermodynam- 
ics in low dimensions by controlling the boson number. 
Constraining the total magnetization to be zero, Taka- 
hashi [^4| obtained an excellent description of the low- 
temperature thermodynamics of one-dimensional Heisen- 
berg ferromagnets. His idea that the thermal spin devi- 
ation, that is, the number of thermally induced bosons, 
should be equal to the ground-state magnetization may 
be replaced by (21j 
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FIG. 2. Dispersion relations of the ferromagnetic (o, o) and 
antiferromagnetic (x) spin waves, which lie in the subspace 
of M = 9 and that of M = 11, respectively, for the parameter 
sets (a) and (b). 
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J2 \^(k)\ 2 = 8S + 4s, 



(4) 



i=l,2,3 



for the ferrimagnetic Mnl2 cluster, where n k j is ex- 
pressed as J27 1 ,n 2 ,n 3 =o n 3 P k(ni,n 2 ,n 3 ) with P k being 
the probability of rij spin waves of mode- j appear- 
ing in the fc-momentum state. Equation (0) claims 
that the thermal fluctuation should cancel the staggered 
magnetization instead of the uniform one, in response 
to the present ferrimagnetic ground state. Minimiz- 
ing the free energy F = E g + J2k X^=i,2,3 ^k,j^j(k) + 
fc BTEfcE ni ,n 2 ,„3 Pk(ni,n 2 ,n 3 )\nP k (n 1 ,n 2 ,n 3 ) with re- 
spect to P k at each k under the condition (^) together 
with the trivial constraints S nijT , 2i7l3 Pk{n\,n 2 ,n 3 ) = 1, 
we obtain the optimum distribution functions as 



(5) 



where /i is a Lagrange multiplier due to Eq. (|J). Once 
the distribution is determined, we can readily calcu- 
late any thermal quantities such as the internal energy 
U = E s + J2k 2 3 ™k,jWj(k) and the magnetic sus- 

ceptibility x = [(3Mb) 2 /k B T] Z)j=i,2 3 "-fcjU + n k,j)- 

Although the above-demonstrated modified spin-wave 
scheme of the ferrimagnetic version generally works well 
in low dimensions |2^,|26|], it is not yet applicable to the 
Mnl2 cluster as it is. Due to the significant Jahn- Teller 
distortion, the lowest excited states of M = ±9 are sepa- 
rated from the ground states of M = ±10 by a finite en- 
ergy A, which is incompatible with the condition (Q). In 
isotropic ferrimagnets, there exists a zero-energy excita- 
tion and therefore a certain number of bosons naturally 
survive at low temperatures. The grandcanonical con- 
straint (|j) not only works so as to suppress the thermal 
divergence of the boson number at high temperatures but 
also gives a precise description of the low-temperature 
thermodynamics ]25( ]. On the other hand, once a gap A 
opens, the boson number should exponentially decreases 
as ex e~ A / kl>T at low temperatures, while the constraint 
(|]) still keeps it finite even at T — ► 0. In order to elimi- 
nate the shortcoming, we replace Eq. fl) by 



E E R *J E hMfc)| 2 = (8S + 4 S K 



A/k B T 



(6) 
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FIG. 3. Temperature dependences of the zero-field mag- 
netic susceptibility for the parameter sets (a) (dotted line) 
and (b) (solid line). Experimental observations [5,9] (o) are 
also shown for reference, where we have assumed that g — 2 
[13]. 

This is quite natural modification of the theory, because 
the new constraint @ remains the same as the autho- 
rized one (|j) except for the sufficiently low-temperature 
region k^T < A. It is also convincing that Eq. (j^) 
smoothly turns into Eq. (||) as A — > 0. In fact, the thus- 
calculated susceptibility looks reasonable in every aspect. 
In Fig. H, we make a logarithmic plot of the zero-field sus- 
ceptibility as a function of T in order to elucidate its low- 
temperature behavior. Regardless of parametrization, x 
exhibits an initial exponential behavior oc e~ A / k ^ T with 
A ~ 13.7 K. The calculations are further consistent with 
the experimental findings ||[| , implying that the param- 
eter set (b) may better describe the Mnl2 cluster. 
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FIG. 4. Semilog plots of l/T^ as a function of T under no 
external field for the parameter sets (a) and (b). The calcula- 
tions, dashed [Mn(l)], dotted [Mn(2)j, and solid [Mn(3)] lines, 
are compared with experiments [19], o [Mn(l)], o [Mn(2)], and 
x [Mn(3)]. 




FIG. 5. Semilog plots of as a function of H at 

T — 1.4K for the parameter sets (a) and (b). The calcula- 
tions, dashed [Mn(l)], dotted [Mn(2)], and solid [Mn(3)] lines, 
are compared with experiments [19], [Mn(l)], o [Mn(2)], and 
x [Mn(3)]. 



Considering the electronic-nuclear energy-conservation 
requirement, the Raman process should play a leading 
role in the nuclear spin-lattice relaxation. The large zero- 
field energy splitting, which is associated with the single- 
ion anisotropy, not only makes the direct process irrele- 
vant but also reduces the possibility of an odd-number- 
of-magnon process being realized. The Raman relaxation 
rate for the 55 Mn(i) nucleus is given by 



1 = 47Tft(ff/X B 7N) 2 p -E n /k B T 

x I (mlAiaf An)] 5(E m - E n - Tk^ 



(7) 



where erf i = sf,Si,Sf for i — 1,2,3, respectively, A t is 



the dipolar coupling constant between the nuclear and 
electronic spins on the Mn(i) site, wn = 7N-ff is the Lar- 
mor frequency of the nucleus with 7n being the gyro- 
magnetic ratio, and the summation is taken over all 
the electronic eigenstates \n) with energy E n . Taking ac- 
count of the significant difference between the electronic 
and nuclear energy scales (fiui-^ < 10~ 5 J), the relaxation 
rate (]?]) is expressed in terms of the spin waves as 



1 



T, 



(i) 



N 



J2 E 



1 j = l,2,3{fc,fe'} 

x \i>ij (fc) 1 2 (k') \ 2 n k j (fife/ j - 



(8) 



where X]{fe k'} denotes the limited summation of (k,k') 
over (-tt.-tt), (-§,-§), (0,0), (§,§), (§,-§), and 
(-§,§), and Aojj(k) = [uj(k + 2%/N) - Wj(k)]/(2n/N). 

We show the thus-calculated l/T^ as functions of T 
(Fig. ^) and H (Fig. (bJ) in comparison with the measure- 
ments [jl9j , where the coupling constants A4 are the only 
adjustable parameters in reproducing both the T- and In- 
dependences and have been chosen as Table Q. Though 
the calculations (a) and (b) look similar, Table | clearly 
shows that the parameter assignment (b) much more rea- 
sonably describes the Mnl2 cluster. In response to ask- 
ing whether the exchange interaction J3 is ferromagnetic 
or antiferromagnetic |jl]], we definitely answer 
that it is antiferromagnetic. Furthermore, we are skepti- 
cal of neglecting J4 |],[l(|, for which the S — 10 ground 
state is much less stable [[14] and the coupling constants 
are significantly underestimated in our calculation. 

The observed temperature dependence of 1 /t[ 1%> , which 
is the same for i — 1,2,3, suggests that the spin dynam- 
ics of the local Mn moments is completely correlated. All 
the observations can indeed be interpreted within the Ra- 
man relaxation process. In the momentum summation 
in Eq. (||), the contribution from (j;k, k') = (1;0, 0) 
(S = 10, M — 9) is predominant, while those from 
(j;|fc|,|fc'|) = (1, §,§) {S = 9,M = 9) at most amount 
to a few percent of the total at T — 1.4 K. The spin 
dynamics is thus confined within fluctuations of the to- 
tal spin moment S = 10 of the ground state at suffi- 
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ciently low temperatures. Although such a picture is in- 
tuitively convincing, we have to pay attention to a recent 
experimental report (27) that the hyperfine field of the 
Mn 3+ ion is in fact anisotropic and exhibits a predomi- 
nant dipolar contribution, whereas that of the Mn 4+ ion 
is isotropic and originates from the Fermi contact. In 
this context, it is interesting to compare carefully the 
theoretical (Af 1 ) and experimental (Af x ) findings for the 
coupling constants. Assuming the set (b), Af 1 ~ 2.4Af x , 
Af ~ 1.7 Af, and Af ~ 1.2Af. Somewhat larger de- 
viation of the theory from the experiment for A\ implies 
that the nuclear spin-lattice relaxation on the Mn(l) site 
may not primarily be Raman active but be strongly in- 
fluenced by the surrounding Mn 3+ ions. 

The slight difference between the field dependences 

(2) (3) 

of l/Tp and 1/Tp can not be elucidated within the 
present calculation, but this is probably due to the equal 
treatment of the Mn(2) and Mn(3) sites, <?2 = <?3 an d 
D2 = D 3 . Considering that experimental analyses [ p~8|Jl9) 
have not yet entered into such details, our first step to- 
ward the microscopic description of the Mnl2 cluster is 
really successful. We stress that the success of our spin- 
wave description contributes in itself toward verifying the 
recently reported intramolecular spin derealization p5[ . 
Our theory promises future investigations into various 
molecular magnets over the static, dynamic, quantum, 
and thermal properties. 
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TABLE I. Estimates of (g/iBTi-Mi) 2 m the unit of 
(rad • Hz) 2 under the parametrizations (a) and (b) compared 
with experimental findings [19]. 





Mn(l) 


Mn(2) 


Mn(3) 


Parameters (a) 


7.9 x 10 17 


3.2 x 10 18 


3.7 x 10 18 


Parameters (b) 


6.5 x 10 16 


2.7 x 10 17 


3.0 x 10 17 


Experimental 


1.1 x 10 16 


9.3 x 10 16 


2.0 x 10 17 
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